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Abstract 

The variation of radii, lifetimes, correlation lengths and the chaoticity para- 
meter of particle sources with the collision energy -^s is studied by analysing and 
comparing two-particle Bose-Einstein correlation data obtained by the NA22 
Collaboration at ^/s = 22 GeV and by the UAl-Minimum-Bias Collaboration 

0^ at ^/s = 630 GeV . The UAl-data are found to be inconsistent with a static 

source, whereas both the data at ^/s = 22 GeV and at 630 GeV can be descri- 

^ bed in an expanding source model. In the latter model, the s-dependence of the 

data implies an increase of the transverse radius with s and a decrease of the 

correlation lengths in longitudinal and in transverse direction. The chaoticity 

^^ parameter increases or remains approximately constant as a function of the CM 

energy. 






O 



^r^ The study of the energy dependence of physical quantities has been an important 

■^ tool in multiparticle dynamics. Thus e.g. the discovery of the violation of KNO 

,__! scaling with the increase of energy from ISR to SPS energies has been one of the 

■^ most stimulating findings in soft strong interactions in the last decade. 

While it is now understood that this phenomenon is due to the energy dependence 
of long range correlations [1] it was unclear whether short range correlations (SRC) do 
also depend on energy. Besides its heuristic interest this question is also of practical 



o^ 



Oh 



Q interest since it is related to the energy dependence of radii and lifetimes which 

1-^ determine via Bose-Einstein correlations (BEC) the SRC of identical particles. 

Here we present the first results of an investigation of the energy dependence of 
source parameters in BEC. These results were made possible by 

• the publication of BEC data at two quite different values of ^/s, i.e. ^/s = 
22GeV [2] and ^ = 630GeV [3] 
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• the development of a space-time theory of BEC [4]. 

Two models were considered for this analysis (for details on these models cf. [5]). 

1. A static source defined by the source distributions fch (for the chaotic part), fc 
(for the coherent part) and the correlation C: 
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where Rch,a and Rc^a{oi = 0,_L, ||) are the lifetimes, transverse radii and longi- 
tudinal radii of the chaotic source and of the coherent source, respectively, and 
where La{a = 0,_L, ||) are the correlation time and the corresponding correla- 
tion lengths in transverse and in longitudinal direction. To specify the relative 
contributions of the chaotic and the coherent component, one needs to fix the 
value of the (momentum dependent) chaoticity parameter p[k^) at some arbi- 
trary scale, e.g., po = p[k^ = 0). 



2. A longitudinally expanding source 
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T is the longitudinal proper time and rj the space-time rapidity. 

The number of free parameters can be reduced by relating them to the shape of 
the one-particle inclusive distribution pi{k). For the static source, one has 
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Figure 1: Best fits to the UAl data 



and for the expanding source, 
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where iV is a normahzation factor and 
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For the static source, a plateau in rapidity requires i?c,o = Rc,\\ = Lq = Z|| = 0. 
Furthermore we make the simphfying assumption that the particle spectra have the 
same shape for the chaotic and for the coherent component, i.e., Rl± = Rc,±- This 
implies that the fraction of chaotically produced particles is momentum independent. 
Therefore, in what follows we have dropped the subscript "0" in the notation for 
the chaoticity parameter. Finally, the value of Rl± = Rc,± can be determined by 
requiring the average transverse momentum {k±) as calculated from Eqs. (7,8) to 
agree with the experimentally observed value^. 

Thus, the remaining four free parameters are i?c/j,o, Rch,\\7 Rch,± and p for the 
static source, and Tq, Rch,±7 L^ and p for the expanding source. 



^A comparison of NA22- [6] and Tevatron-data [7] shows that the average transverse momentum 
of pions depends only weakly on the CM energy. In our calculations, we have taken {k±) ~ 0.37 GeV . 
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Figure 2: Curves for the NA22 phase space cuts, obtained for the same parameter 
values as in Fig. la. 



Since the data are expressed in terms of the invariant variable Q^, one has to 
calculate the correlation function in terms of this variable. The correlation function 
reads^: 



C2(g' 
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Below we discuss our results: 

Fig. 1 shows the UAl-data"^ for the two-particle BEC as a function of Q^ together 

with the fits obtained in refs. [8, 9] for a purely chaotic source [p = 1), for the static 

and the expanding source, respectively. The corresponding parameter values and the 



^It is worth mentioning that the same type of integrals occur in the interpretation of "intermitt- 
ency" data. In ref. [8, 9] it was shown that a conventional source of the type defined above with 
a fixed geometry (at a given energy) can explain all the so-called intermittency observations, i.e. 
the dependence of moments on Q and there is no need for a fractal source as suggested in [10]. As 
a matter of fact a fractal source does not imply even a simpler explanation of the data because it 
involves an exotic length scale of the order of 6-7 fm for which there is no simple interpretation 

^We have normalized both the UAl-data [3] and the NA22-data [2] to unity at Q^ = 1 GeV"^ in 
order to subtract the effects of long range correlations which were shown to play an important role 
at UAl-energies[l]. 



Static source 


V^ [GeV] 


P 


Rch,o [fm] 


Rch,\\ [fm] 


Rch,± [fm] 


xyd.o.f. 


630 


1.00 


1.00 


3.50 


1.05 


1.65 


22 


1.00 


3.80 
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1.00 


1.42 


22 
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1.00 


1.00 



Table 1: Fit parameter values and x^/d-o.f. for the static source. 



Expanding source 


V^ [GeV] 


P 


To f m 


Rch,± [fm] 


L,^ [fm] 


xyd.o.f. 


630 


1.00 


2.00 


1.00 


0.40 


1.33 


22 


1.00 


1.90 


0.65 


0.90 


0.99 


22 


0.60 


1.90 


0.70 


0.90 


0.90 



Table 2: Fit parameter values and x^ I d-o.f . for the expanding source. 



values of x^ /d.o.f. are listed in Tables 1 and 2. 

The NA22-data are plotted in Fig. 2. Clearly, the BEC function C2{Q^) at 
^/s = 22 GeV is broader than the one at 630 GeV. Naively, one might conclude 
from this observation that the radii and/or lifetimes of the source increase with CM 
energy. However, one must bear in mind that the correlation functions were obtained 
in different phase space regions which enter into the integrations in (10): for the 
NA22-data, is defined through the cuts \y\ < 2^ k± > 0, whereas for the UAl-data 
one has \y\ < 3^ k± > 0.15 GeV. We have applied the fit parameters used for the 
description of the UAl-data (Fig. 1) to calculate C2{Q^) in the kinematical region 
that corresponds to the NA22-data. The results, both for the static source and for 
the expanding source, are shown in Fig. 2. Obviously, they fail to reproduce the 
NA22-data. This implies that the observed difference in the widths of C2{Q^) in 
the NA22- and in the UAl-experiment cannot be explained through the difference 
in the rapidity and transverse momentum cuts alone. Rather, it refiects a genuine 
dependence of the BEC on the CM energy ^/s of the collision. 

To describe the NA22-data, one needs to choose a set of parameter values different 
from the ones used to fit the UAl-data. The resultant space-time parameters are listed 
in Tables 1 and 2, and the corresponding correlation functions are compared to the 
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Figure 3: Best fits to the NA22 data 



NA22-data in Fig. 3. 

Up to this point, we have always assumed that the particles are emitted from 
a purely chaotic source. Before we can discuss the quality of the fits and the s- 
dependence of the quantum statistical parameters of the source, we have to address 
the question if the data allow any conclusions concerning the presence of a coherent 
component and its energy dependence. 

The x^ /d.o.f. for the best fits are plotted as a function of the chaoticity p in Figs. 
4a and 4b, for the UAl-data and the NA22-data, respectively, where the dashed 
curves correspond to the static source and the solid curves to the expanding source. 
According to the x^-criterion, parameter sets that correspond to points above the 
straight line labeled 5% (1%) can be ruled out at the 5% (1%) level, i.e., the probabi- 
lity that the data are consistent with the model is below 5% (!%)• At ^/s = 630 GeV , 
the static source can be ruled out on the 5% level. On the other hand, the expanding 
source can fit the data with reasonable values of x^ /d.o.f. for chaoticities in the range 
0.9 < p < 1. For the data at ^/s = 22 GeV , both the static and the expanding source 
models are acceptable for a wide range of chaoticity parameters (0.37 < p < 1 for 
the static source and 0.31 < p < 1 for the expanding source). As an example, we 
have listed the values of the fit parameters and x^ /d-o.f. for a chaoticity of 0.6 at 
^/s = 22 GeV in Tables 1 and 2. Note, however, that even at NA22-energies one 
might have some doubts about the static source model as the fits require lifetimes of 
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Figure 4: x^/c/.o./. for the best fits as a function of the chaoticity parameter p, (a) 
for the UAl-data and (b) for the NA22-data. As before, the dashed curves are for 
the static and the sohd curves for the expanding source modeL Points above the 
horizontal hne labeled 5% (1%) correspond to chaoticity parameter values with a 
probability below 5% (1%) for the data to be consistent with the model 



3 to 4 fm^ i.e., values which are suspiciously large for hadronic collisions. Neverthe- 
less, in principle one cannot exclude the possibility that particles are emitted from a 
static source at ^/s = 22 GeV and from an expanding source at ^/s = 630 GeV. As 
was pointed out in [9], a sensitive test which would allow a clear distinction between 
these two types of sources is to measure C2 as a function of the two variables Qq and 
q^. There are no published data in these two variables at ^/s = 22 GeV and at 
^/s = 630 GeV , and experimentalists are urged to fill this gap. 

Although a static source at y/s = 22 GeV cannot be completely ruled out at this 
point, in what follows we shall restrict our discussion of the energy dependence of the 
quantum statistical parameters to the expanding source model. As we have seen, the 
data at ^/s = 22 GeV are consistent with chaoticity parameters in the range between 
0.3 and 1. Fortunately, for the expanding source at 22 GeV the values of the other 
fit parameters are not very sensitive to the value of p. As can be seen in Table 2, for 
p = 1 one has Rch,i. = 0.65 fm, L^ = 0.90 fm and Tq = 1.90 fm^ while for p = 0.6 the 
values are Rch,± = 0.70 fm^ L^ = 0.90 fm and Tq = 1.90 fm. This enables us to draw 
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Figure 5: Predictions for C2{Q^) in the momentum space region \y\ < 2, 
k± > 0.15 GeV , for the expanding source at ^/s = 22 GeV and at ^/s = 630 GeV . 
The fit parameter values are the same as in Figs. 1 and 3. 



conclusions about the s-dependence of i?c/j,_L, Ln and Tq although the NA22-data do 
not yield strong constraints on the value of p. Going from NA22- to UAl-energies, the 
longitudinal proper time Tq is found to be almost unaffected (it changes from 1.9 fm to 
2 fm). On the other hand, the transverse radius Rch,i. increases by about 40% (from 
~ 0.7 fm to 1 fmi) and the correlation length in longitudinal direction Z^ is reduced 
by more than 50% (it drops from 0.9 fmi to 0.4 fmi). Concerning the chaoticity 
parameter one can conclude that it either increases with s or remains approximately'* 
constant. The correlation length in transverse direction L_\_ decreases from the value 
of 0.8 fmi at 22 GeV to 0.6 fmi at 630 GeV . If there is a coherent component, the 
transverse radius of the coherent source remains constant as a function of s, at a value 
of i?c,_L ~ 0.5 fmi. 

Of course it would be interesting to compare the BEC data at ^/s = 22 GeV and 
at ^/s = 630 GeV in the same momentum space region. The kinematical cuts are 
\y\ < 2, A;_L > for the NA22-data and \y\ < 3, A;_l > 0.15 GeV for the UAl-data. We 
would therefore urge both experimental groups to re-analyze their data and present 
C2{Q^) in the restricted phase space region \y\ < 2^ k± > 0.15 GeV . This would allow 



■*The data are also consistent with a slight decrease from p 
/i= 630 GeV, cf. Fig. 4. 



1 at yi = 22 GeV to p ~ 0.9 at 



to explicitly see the s-dependence of the correlation function without having to take 
into account the effects of different phase space cuts. Fig. 5 shows our predictions 
for C2{Q^) in the region \y\ < 2^ k± > 0.15 GeV for the expanding source model, at 
^/s = 22 GeV and at ^/s = 630 GeV. The parameter values of Tq, Rch,±7 L^ and 
p = 1 correspond to the fits to the data in Figs. 1 and 3. As can be seen by comparing 
Fig. 5 with Fig. 1, for the expanding source the restriction in rapidity (from |j/| < 3 
to \y\ < 2) leaves the shape of the correlation function practically unchanged (cf. 
also [9]). On the other hand, a comparison of Fig. 3 and Fig. 5 shows that the 
cut in transverse momentum (from A;_l > to A;_l > 0.15 GeV) leads to a significant 
increase of the width of the correlation function. According to our prediction, the 
BEC function C2{Q^) at ^/s = 22 GeV will be considerably broader than the one 
at ^/s = 630 GeV ^ if the data at both energies are taken in the same momentum 
space region. This refiects the increase of the transverse radius and the decrease of 
the correlation lengths as a function of s. 

Note that here we have not treated separately the contributions from the different 
resonance decays^. Thus the radii and lifetimes extracted from the BEC data and 
listed in Tables 1 and 2 must be regarded as effective radii and lifetimes in so far as 
they represent an average of the direct pion component and the components related 
to pions produced from resonance decays. We have also not attempted to correct for 
the effects of particle misidentification. 

Our results concerning the s-dependence of the space-time parameters of the 
expanding source are in qualitative agreement with the experimental observation 
[12, 13, 14, 15] that at fixed ^/s the width of the BEC function decreases with incre- 
asing rapidity density dn/dy. Since the average dn/dy increases with s, the dn/dy- 
dependence of the BEC data suggests that the width of the correlation function 
decreases, which in the expanding source model is refiected in an increase of Rch,± 
and a decrease of Z^, L± with increasing CM energy. Unfortunately, the "effective 
radii" extracted from the data at fixed dn/dy do not distinguish between components 
in beam direction and in the transverse plane. While the s-dependence of the space- 
time parameters agrees with what one would expect from the behaviour of the data 
as a function of dn/dy ^ this is not the case for the chaoticity parameter p. The data 
suggest a decrease of the A-parameter with increasing dn/dy. While A does not agree 
with the quantum statistical chaoticity p), qualitatively a decrease of A translates 
into a decrease of p. On the other hand, in our present analysis of the data we have 
found that p either increases or remains approximately constant as the CM energy 
increases from NA22- to UAl-energies. The resolution of this apparent discrepancy 
may be related to the fact that the extraction of a reliable value for the chaoticity 
parameter requires knowledge of the behaviour of the BEC function near the origin, 
i.e., at small momentum differences. It would therefore be useful also from this point 
of view if the experimentalists could measure C2{Q^) at fixed dn/dy down to values 



®A formalism which explicitly takes into account the effect of resonance decays on BEC is pre- 
sented in ref. [11]. 



of Q — 30 MeV. The analysis of such data could help to determine the dependence 
of the quantum statistical parameters on the two variables dn/dy and s. 

In ref. [16], the dependence of the normalized factorial moments of the multipli- 
city distribution on the width of the rapidity interval was analysed for NA22- and 
for UAl-data. It was found that both the chaoticity parameter and the correlation 
length in rapidity ^ increase with s. As an increase of ^ corresponds to an increase 
of Lrj in the present formalism, the s-dependence of ^ obtained in [16] would appear 
to be in contradiction with the s-dependence of L^ obtained here. However, a note 
of caution is in order. Our present results cannot be compared to those of ref. [16] 
in a straightforward manner for a number of reasons. Firstly, the calculations in [16] 
were based on an ansatz for the correlation function in one variable only (rapidity). 
Secondly, at that time the factorial moments at Collider energies were known only 
for charged particles. Finally, and most importantly, in Ref. [16] it was assumed that 
only short-range correlations (SRC) contribute to the correlation function. However, 
recently it has been found [1] that at Collider energies there is a significant contri- 
bution due to long-range correlations (LRC), and that LRC increase strongly with 
s. Indeed it appears that in ref. [16] the strong increase of fiuctuations due to the 
increase of LRC was attributed to an increase of SRC and hence, an increase of ^ 
with the CM energy. 

On the theoretical side, there exist surprisingly few predictions concerning the 
dependence of BEC on the CM energy ^/s. While the Lund model [17] apparently has 
no s dependence for BEC a decrease of the width of the BEC function in longitudinal 
direction, 6*2(511), with increasing s was predicted by Bowler in the context of a string 
model[18]. A similar behaviour is also expected in a hydrodynamic scenario: for 
a longitudinally expanding source, the width of 6*2(511) i^ approximately inversely 
proportional to the freeze-out proper time Tf [19]. If the formation time Tq of the 
thermalized matter is assumed to be s-independent, Tf increases with the initial 
energy density which in turn increases with s. Consequently, the width of 6*2(511) 
decreases with increasing CM energy. 

An increase of the transverse radius with s can be understood in the context 
of a simple geometric description of multiparticle production. From ^/s = 22 GeV 
to 630 GeV the total pp (pp) cross section atot increases by about 60% (cf. the 
parametrization of crtot{-s) in ref. [20]). If one assumes that the transverse source 
radius scales approximately like a^^^ one would predict an increase of Rj_ by ~ 30%. 
This is to be compared with the value of ~ 40% found in our analysis based on 
the model of a longitudinally expanding source^. For a hydrodynamically expanding 
source the s-dependence of the longitudinal and transverse fiow components may also 



^Indeed, the fact that an increase of R± is to be expected from simple phenomenological consi- 
derations can be regarded as further support for the expanding source model and evidence against 
the validity of the static source model. For the case of a static source at ^/s = 22 GeV, our analysis 
of the data would imply that the transverse radius takes approximately the same value at 22 GeV 
and at 630 GeV (cf. Tables 1 and 2). 



contribute to the increase of the transverse radius. As was shown in ref. [21], for 
reasonable values of the freeze-out temperature transverse flow can lead to a decrease 
of the transverse source radius extracted from BEC data (cf. Fig. 3 in [21]). At 
higher energies, the longitudinal expansion dominates and transverse flow becomes 
less important (cf. ref. [22]). Therefore, the transverse source radius can be expected 
to increase with the CM energy. 

The fact that the observed change with energy of source geometry is in agreement 
with expectations derived from independent considerations proves that the theory of 
BEC with its space-time formulation [4] has come out of age and gives confldence for 
its use in applications like quark-gluon plasma search. This conclusion is reinforced by 
the good agreement of "intermittency" data with results obtained from Bose- Einstein 
correlations [8, 9]. 

This work was supported by the Federal Minister of Research and Technology 
under contract 06MR731 and the Deutsche Forschungsgemeinschaft. We are indebted 
to B. Buschbeck (UAl Collaboration), W. Kittel (NA22 Collaboration) and U. Ornik 
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